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Abstract 

The  goal  is  to  fulfill  the  US  Navy’s  need  for  a  relocatable,  robust  operational 
coastal  forecast  system  by  developing  and  transitioning  a  high-resolution,  coastal 
circulation  model  into  the  Naval  Oceanographic  Office’s  operational  environment. 
The  coastal  circulation  model  undergoing  transition  is  a  three-dimensional, 
finite-element  based  hydrodynamic  model,  the  Advanced  Circulation  Model  for 
Shelves,  Coastal  Seas,  and  Estuaries  (ADCIRC).  Its  unstructured  grid  allows 
modeling  complex  coastal  regions  at  fine  spatial  scale.  Scripts  and  programs  arc 
developed  to  automate  major  tasks  of  the  forecast  run  stream  such  as  initial  setup, 
forcing  data  acquisition,  model  configuration  and  post-processing.  An  application 
of  the  ADCIRC-based  forecast  system  to  the  Korean  West  Coast  region  is 
presented. 


I.  INTRODUCTION 

Navy  activities  are  global,  typically  in  non-US  and  denied  waters;  operators  have 
limited  time  to  set-up  and  configure  models  and  make  predictions  to  support  missions. 
It  is  critical  to  develop  a  modeling  system  that  can  be  configured  quickly,  and  then 
perform  model  simulations  and  predictions  automatically  with  minimal  human 
intervention.  The  goal  is  to  fulfill  the  US  Navy’s  need  for  a  relocatable,  robust 
operational  coastal  modeling  system  by  developing  and  transitioning  a  high-resolution, 
coastal  circulation  model  into  the  Naval  Oceanographic  Office’s  (NAVO)  operational 
environment. 

A  relocatable  system  simply  means  a  system  can  be  moved  to  different  geographic 
regions  or  scales  easily,  with  minimal  changes  in  domain  configuration.  NRL  has 
developed  relocatable  systems  for  ocean  models  (RELO-NCOM)  [1]  and  tide  models 
(PCTides)  [2]  in  the  past  few  years.  As  the  Navy  increases  operations  in  coastal, 
estuarine  and  shallow  regions,  a  relocatable  coastal  modeling  system  is  needed 
(RELO-ADC1RC)  that  has  the  capability  to  make  accurate  water  level  and  current 
predictions  in  complex  coastal  regions.  Two  major  challenges  in  developing  such 
system  are  open  boundary  specification  and  grid  generation  for  models  based  on 
unstructured  grids.  The  grids  for  several  strategic  geographic  regions  are  generated  in 
advance  and  stored  in  a  grid  database.  An  automated  grid  generation  tool,  MeshGUl  [3], 
is  also  supported  for  areas  not  yet  available  in  the  database.  For  the  surface  forcing, 
meteorological  wind  and  atmospheric  pressure  fields  are  provided  by  the  Navy  Coupled 
Ocean  Atmospheric  Mesoscale  Prediction  System  (COAMPS)  [4]  Scripting  programs 
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are  used  to  integrate  individual  modeling  components  and  to  automate  tasks  such  as 
creating  and  configunng  model  input  files,  acquiring  surface  atmospheric  forcing  data 
and  open  boundary  tidal  information,  checking  model  run  status,  reporting, 
post-processing  and  mapping  final  products.  This  article  describes  the  model 
configuration,  grid  generation  tools,  automated  scripting,  and  an  example  application  of 
the  system. 


II.  ADCIRC  MODEL 

The  core  coastal  circulation  model  used  in  the  forecast  modeling  system  is  a 
three-dimensional,  high-resolution,  finite-element  based  hydrodynamic  model, 
ADvanced  CIRCulation  Model  (ADCIRC).  Its  unstructured  grid  system  allows 
modeling  complex  coastal  regions  at  fine  spatial  scale.  With  properly  specified 
topography/bathymetry,  meteorological  forcing,  tidal  potential  and  boundary  forcing, 
this  model  can  accurately  simulate  water  levels,  currents  and  area  of  inundation  due  to 
wind,  tides  or  other  forcing. 

Major  improvements  have  been  made  since  the  ADCIRC  model  first  developed  in  the 
early  nineties.  The  once  separated  2D  and  3D  versions  of  the  code,  arc  now  contained 
within  a  single  unified  program.  The  code  has  been  translated  into  the  FORTRAN  90 
language  standard  and  is  completely  parallelized  using  the  MP1  protocol  and  metis 
domain  decomposition  software.  The  code  structure  itself  is  modular  with  individual 
subroutines  containing  global  dimensioning,  cold  and  hotstart  initialization,  parameter 
specification,  wind  forcing,  harmonic  analysis,  three-dimensional  calculations,  and  the 
main  time  stepping  loop.  Furthermore,  the  code  has  become  multi-algorithmic  with 
options  for  selecting  various  implementations  of  the  conservative  or  non-conservative 
Generalized  Wave  Continuity  Equation  (GWCE)  and/or  momentum  equation 
formulations.  One  notable  addition  to  the  ADCIRC  v45.I  I  model  code  is  an  improved 
wetting  and  drying  algorithm.  This  implemented  wetting/drying  capability  has 
significantly  improved  realism,  performance  and  stability  of  the  model  when 
wetting/drying  and  inundation  are  important. 

A  complete  description  of  the  theoretical  basis  for  ADCIRC  is  given  by  [5,  6]  which  is 
found  online  with  the  User’s  manual  at 

http://www.adcirc.org/documcnt/ADCl  RC  titlc_pagc.html.  A  comprehensive  validation 
of  the  2D  and  3D  model  response  under  meteorological  forcing,  tides  and  inundation  for 
various  geographic  regions  can  be  found  in  the  newly  released  NRL  Validation  Test 
Report  [7]. 


III.  GRID  DATABASE 

The  system  uses  a  database  concept  to  make  available  pre-existing  unstructured  meshes 
for  the  ADCIRC  numerical  model.  One  of  the  main  features  of  the  ADCIRC  model  is  its 
ability  to  apply  spatially  varying,  high  resolution  finite  element  grids  to  represent 
complex  coastal  regions.  As  a  result,  additional  effort  is  needed  to  create,  edit  and  refine 
the  triangular  meshes  when  compared  with  traditional  finite  difference-based  gridding 
systems.  In  the  database  approach  used  for  this  system,  grids  for  several  strategic 
regions  for  the  Navy,  including  the  West  Korean  Coast  have  been  created,  and  stored  in 
the  database.  An  NRL  in-house  developed  mesh  generation  and  editing  tool  with 
Graphical  User  Interface  (GUI),  MeshGUI,  is  included  in  the  system  in  case  users 
require  the  creation  of  meshes  for  areas  not  yet  available  in  the  existing  grid  database. 
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IV.  SYSTEM  SOFTWARE  DEVELOPMENT 

One  of  the  major  requirements  in  developing  an  operational  modeling  system  is  to 
automate  routine  tasks  as  much  as  possible  while  still  maintaining  system  flexibility  and 
allowing  manual  override.  To  meet  this  requirement,  a  software  program  is  needed  to 
integrate  individual  modeling  components  and  to  automate  tasks  such  as  configuring 
model  input  files,  acquiring  surface  wind  forcing  data  and  open  boundary  tidal 
information,  checking  model  run  status,  reporting,  postprocessing  and  mapping  final 
products.  The  software  developed  in  support  of  system  automation  includes  a 
collection  of  shell  scripts,  awk  scripts,  Perl  programs,  C  programs,  Fortran  programs 
and  third-party  software.  Each  program  is  designed  to  serve  a  specific  purpose  or 
accomplish  a  particular  task.  Figure  1  shows  the  schematic  diagram  of  this  modular 
system 

All  the  programs  and  scripts  are  grouped  into  five  modules  according  to  their  functions: 
1)  System  Setup,  2)  Data  Acquisition,  3)  Model  Configuration,  4)  Model  Run  and  5) 
Products  Generation  and  Dissemination.  All  of  the  scripts  and  modules  are  designed  to 
be  manually  invoked  from  the  command  line  if  errors  are  found  or  if  a  scheduled  task  is 
not  completed.  This  feature  also  allows  manual  processing  to  be  performed  for  cases 
where  input  files  at  the  time  of  execution  are  missing  or  incomplete.  Main  functions  and 
programs  in  each  module  are  described  as  follows: 

Module  1:  System  Setup 

This  module  contains  programs  and  scripts  that  setup  the  operating  system’s 
environmental  variables,  create  proper  file  directory  structures,  and  initialize  starting 
and  ending  timestamps. 

Module  2:  Data  Acquisition 

All  the  required  data  for  the  modeling  system  are  collected  and  processed  in  this  module. 
Software  used  for  downloading  data  from  remote  servers  and  reformatting  data  into 
ADCIRC-compatible  files  for  model  runs,  are  placed  in  this  module.  Required  input 
files  for  the  ADCIRC  modeling  system  include  meteorological  wind  forcing  (fort.22), 
open  boundary  tidal  information  (fort.  15)  and  river  inflows  (fort. 20)  (currently  not 
used).  As  an  example,  the  u,  v  components  of  wind  and  pressure  fields  are  extracted 
from  the  Western  Pacific  region  (WPAC)  COAMPS  standard  27-km  data  at  3-hour 
intervals.  An  option  for  accepting  nested  high  resolution  COAMPS  fields  (9km  or  3km) 
will  be  provided  Two  interpolation  programs  are  then  applied  to  create  spatially  and 
temporally  appropriate  and  continuous  high  resolution  meteorological  fields. 

Although  not  currently  used,  a  script  is  written  to  allow  manual  extraction  of  surface 
elevation  and  current  components  from  the  Navy  Coastal  Ocean  Model  (NCOM)  [8] 
model  outputs.  Once  the  unstructured  mesh  is  finalized,  locations  of  mesh  nodes  along 
the  seaward  boundary  can  be  extracted  and  saved  in  an  ASCII  text  file  that  is  then  used 
to  identify  exact  locations  within  the  NCOM  grid  to  extract  the  NCOM  parameters  of 
interest  for  the  analysis  time  and  forecasts  desired.  This  can  be  used  to  provide  surface 
elevation  or  normal  flux  boundary  conditions  for  ADCIRC  forcing. 
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Module  3:  Model  Configuration 

This  module  consists  of  programs  and  shell  scripts  that  create  necessary  ADC1RC 
model  related  input  files,  fort.  14,  fort.  15,  fort.20  and  fort.22.  In  addition,  scripts  that  run 
the  model  preparation,  execution  and  post-processing  code,  adcprep ,  padcirc  and 
adcpost ,  respectively,  are  placed  in  this  module.  The  four  required  input  files  include  the 
fort.  14  grid  file,  the  fort.  1 5  model  parameter  file,  the  fort.20  flux  boundary  forcing  file, 
and  the  fort.22  surface  meteorological  forcing  file  are  prepared  differently  as  described 
below: 

The  grid/bathymetry  file  is  copied  from  the  grid  database  to  the  local  directory  and 
renamed  fort.  14  to  interface  directly  with  ADCIRC  model  input  file  units.  Fort.  15  is 
the  main  model  configuration  and  parameter  file.  Fixed  parameters  are  save  in  a 
template  fort.  15  file  while  run  specific  variables  such  as  starting  time,  run  period  or  user 
specified  parameters  are  extracted  from  Module  1  and  appended  to  the  template  file. 
This  file  needs  to  be  updated  for  each  simulation.  Fort.20  is  the  river  inflow  forcing  but 
is  not  currently  supported  in  this  implementation  of  the  model  forecast  system.  The 
fort.22  file  is  created  based  the  data  interpolated  from  the  aggregated  COAMPS  wind 
data  acquired  from  Module  2. 

Three  additional  scripts  that  setup,  run  and  post-process  parallel  ADCIRC  model  runs 
are  also  included  in  this  module.  The  a deprep  program  prepares  parallel  runs,  performs 
domain  decomposition,  and  creates  necessary  sub-directories  based  on  the  user 
specified  number  of  processors  and  copies  necessary  files  to  proper  directories.  The 
model  executable,  padcirc ,  runs  the  ADCIRC  model  as  a  parallel  program  using  the 
Message  Passing  Interface  (MP1)  software  package  (OpenMPl).  Once  model  execution 
is  complete,  the  post-processing  program  adcpost,  is  executed  to  put  all  regionalized 
output  into  global  files.  Default  input  template  files  are  also  written  to  provide  NAVO 
users  additional  flexibility  to  run  the  model  in  an  automated  batch  mode  or  interactively 
in  the  case  where  the  hardware  platform  or  the  number  of  processors  needs  to  be 
changed. 

Module  4:  Model  Run 

The  reason  for  separating  Model  Run  as  an  individual  module  is  again  to  provide 
NAVO  operators  additional  flexibility.  In  an  operational  modeling  environment,  it  is 
quite  common  that  model  configuration,  data  acquisition  and  model  simulations  are  all 
done  on  different  hardware  platforms,  for  example,  configuring  model  parameter  files 
on  a  PC  desktop,  processing  model  input  data  streams  on  a  Linux  workstation  while 
running  the  code  on  a  supercomputer  or  parallel  clusters. 

Module  5:  Product  Generation  and  Dissemination 

This  module  contains  software  and  scripts  used  to  generate,  disseminate  and  archive 
map  products  and  model  outputs.  The  product  mapping  scripts  invoke  appropriate 
programs  and  mapping  routines  based  on  users’  preferences  of  mapping  tools  such  as 
the  Generic  Mapping  Tools  (GMT),  MATLAB  or  1DL.  Mapped  products  will  then  be 
place  on  a  designated  web  server  for  the  end  users.  The  planned  products  include 
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graphical  maps  of  forecast  wind  fields  (speed  and  directions),  water  elevations  in  color 
contours  and  currents  with  arrows  showing  directions  and  color  indicating  magnitudes. 

Time  series  of  water  levels  and  currents  at  specific  locations  will  also  included  in  the 
products  suite.  Statistical  metrics  such  as  root  mean  square  error  (RSME),  mean 
absolute  error  (MAE)  and  standard  deviation  (SD)  will  be  computed  based  on  the  user’s 
specification. 


V.  APPLICATION:  WEST  KOREAN  COAST 

An  application  of  the  developed  forecast  system  to  the  West  Korean  Coast  is  used  to 
demonstrate  the  automated  and  operational  capabilities  of  the  modeling  system.  The 
West  Korea  Coast  is  a  highly  dynamic,  tidally-driven  region  with  complex  coastlines 
and  tidal  flat  zones;  it  is  therefore  an  excellent  testbed  for  a  high  resolution  system  that 
has  the  capability  to  model  complex  geometry  and  wet/dry  situations. 

The  grid  bathymetry  was  derived  from  the  Naval  Research  Laboratory  Global  Digital 
Bathymetry  Data  Base  (DBDB2)  2-minute  resolution  dataset 
(http://www7320.nrlssc.navy.mil/DBDB2  WWW).  This  dataset  also  incorporates  the 
high  resolution  Yellow  Sea/Korean  Coast  data  from  Choi  [9].  The  model  grid  domain 
encompasses  the  Yellow  Sea,  the  Bohai  Sea,  and  part  of  the  East  China  Sea  where  the 
open  ocean  boundary  follows  the  200  m  isobaths  between  Taiwan  and  Japan  [101. 
Figure  2  shows  the  Yellow  Sea  regional  grid  that  consists  of  26,752  nodes  with  highest 
resolution  of  420  meters.  The  tidal  potential  and  tidal  constituents  are  extracted  from  the 
FEM99  global  tidal  database  [11,  12]  and  the  present  setup  includes  eight  main  tidal 
constituents  (Ql,  OI,  Kl,  N2,  M2,  K2,  S2  and  PI).  The  Korean  Application  has  been 
running  semi-operationally  on  NRL’s  linux  cluster  since  September  1st,  2009.  The 
current  setting  in  the  model  configuration  is  to  make  48-hour  forecasts  valid  at  0Z.  A 
typical  48-hour  forecast  run  cycle  takes  approximately  8  minutes  wall  clock  time  on  a 
16  processor  linux  cluster.  Current  model  products  include  hourly  2-dimensional  water 
level  and  current  plots  for  the  Korean  Coast  region  and  time  series  of  water  levels  at 
specified  locations.  Figures  3  and  4  show  examples  of  the  water  level  and  currents  plots. 

VI.  SUMMARY  AND  FUTURE  ENHANCEMENTS 

The  framew  ork  of  a  relocatable  forecast  system,  RELO-ADC1RC  ,and  an  application  of 
the  West  Korean  Coast  have  been  developed.  The  system  consists  of  numerous  Fortran 
and  C  programs  and  C  shell,  awk  and  Perl  scripts  that  are  categorized  into  five  modules 
according  to  their  task  and  functions.  At  the  time  of  this  writing,  the  West  Korean  Coast 
application  has  been  configured  and  running  semi-operationally  on  NRL’s  linux  cluster. 
It  is  anticipated  that  an  additional  six  months  will  be  used  to  evaluate  system  robustness 
and  validate  model  performance  before  transition  to  NAVO  for  operational  use. 

Future  enhancements  of  the  modeling  system  include: 

1)  Continued  population  the  unstructured  grid  database  for  Navy  geographic 
regions  of  interest. 

2)  Extraction  of  tidal  constituents  at  the  open  ocean  boundary  condition  using 
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the  Oregon  Tide  Inversion  Software  (OTIS)  [13]  . 

3)  Integration  of  boundary  forcing  from  regional  scale  NCOM  domains  for 
improved  nesting/coupling. 

4)  Incorporation  of  a  river  database  to  account  for  the  time-varying  river  influx  to 
improve  predictions  in  estuarine  regions. 

5)  Upgrade  the  ADCIRC  code  to  include  its  fully  3D  baroclinic  capability  with 
NetCDF  format  input/output  support. 

6)  Explore  the  potential  of  deriving  bathymetry/topographic  data  from  a 
remote-sensing  platform  such  as  Shuttle  Radar  Topography  Mission  (SRTM), 
for  rapid  grid/mesh  generation. 
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Figure  1.  Schematic  diagram  of  RELO-ADCIRC  system 


Figure  2.  Yellow  Sea/Korean  Coast  gnd  with  bathymetry 
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Figure  3.  Example  of  water  elevation  plot 
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